Gap Function with Point Nodes in Borocarbide Superconductor YNi 2 B 2 C 
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To determine the superconducting gap function of a borocarbide superconductor YNi2B2C, the c- 
axis thermal conductivity n zz was measured in a magnetic field rotated in various directions relative 
to the crystal axes. The angular variation of n zz in H rotated within the afe-plane shows a peculiar 
fourfold oscillation with narrow cusps. The amplitude of this fourfold oscillation becomes very small 
when H is rotated conically around the c-axis with a tilt angle of 45°. Based on these results, we 
provide the first compelling evidence that the gap function of YNi2B2C has point nodes, which 
are located along the [100] and [010]-directions. This unprecedented gap structure challenges the 
current view on the pairing mechanism and on the unusual superconducting properties of borocarbide 
superconductors. 
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Superconductivity with unconventional pairing sym- 
metry has been a central subject in the physics of su- 
perconductors. The vast majority of superconductors 
have conventional s-wave symmetry with an isotropic 
gap which is independent of directions over the entire 
Fermi surface. In the last two decades unconventional 
superconductivity with different symmetry has been 
found in heavy fermion materials jjj], high-T c cupratesQ, 
ruthenateQ, and organic compounds Q. Unconven- 
tional superconductivity is characterized by anisotropic 
superconducting gap which becomes to zero (nodes) 
along certain crystal directions j5|. Since the nodal 
structure is ultimately related to the pairing interaction, 
its identification is crucial for understanding the pairing 
mechanism. A common feature in the unconventional su- 
perconductors discovered until now is that they all have 
line nodes in the gap functions, which are located paral- 
lel or perpendicular to the basal planes. It is generally 
believed that these nodal structures appear as a result 
of pairing interactions caused by purely electronic ori- 
gin, instead of being mediated by conventional phonons. 
In fact a close relation between the spin fluctuation and 
superconductivity has been discussed in all of these un- 
conventional superconductors. 

It is well known that there can be another type of gap 
singularity; i.e. point nodes in which the gap turns to 
zero at isolated points on the Fermi surface. The point 
node is thought to be present in the ^4-phase of superfluid 
3 He ||. On the other hand, its existence has never been 
reported in any superconductor so far. To our knowledge, 
the one exception is the B-phase of UPt3 in which a pos- 
sibility of the point nodes coexisting with line nodes was 
suggested theoretically , but with no firm experimental 



evidence. Thus it is still an open question whether a su- 
perconductor can have point nodes, and its clarification is 
an important subject for the study of the unconventional 
superconductivity. 

In this Letter, we direct our attention to a new class 
of superconductors; non-magnetic borocarbide supercon- 
ductors LnNi2B2C, Ln=(Y and Lu) ||. These borocar- 
bide superconductors with elevated transition tempera- 
tures are in the clean regime £ <C I, where £ is the co- 
herence length and I the mean free path. The electronic 
structure is essentially 3D and the crystal structure is 
tetragonal. At an early stage, the gap symmetry of these 
materials was considered to be isotropic s-wave, mediated 
by conventional electron-phonon interactions How- 
ever, recent experimental studies, such as specific heat 
[10, [H], thermal conductivity |1^| , Raman scattering 
[ 1 3 1 , and photoemission spectroscopy (l4| on YM2B2C or 
LUNI2B2C have reported a large anisotropic gap function. 
Despite these studies, the detailed structure of the gap 
function remains unknown. Here, we provide a strong 
evidence that the gap function of YNi 2 B 2 C has -point 
nodes, which are located along the [100] and [010] direc- 
tions. These results are based on the angular variation 
of the thermal transport in a magnetic field rotated in 
various directions relative to the crystal axes, which has 
proved to be a powerful probe for determining the low en- 
ergy quasiparticle (QP) excitation including its direction 
|§ |l|, [L7| f[|, |(], ^ H [23|]. This unprecendented 
gap structure challenges the current view on the mecha- 
nism of the superconductivity and on the many unusual 
properties in the superconducting state, both of which 
have been discussed without taking account the presence 
of the point nodes B |1. 
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FIG. 1: (a)T-dependence of the c-axis thermal conductivity 
k zz in zero field. Inset: Log- log plot of the same data be- 
low 1 K. (b)Field dependence of k zz at low temperatures ( 
H\\ [110]). The solid circles represent the data measured by 
sweeping H after zero field cooling and the open circles rep- 
resent the data measured under the field cooling conditions 
at each temperature. 



Single crystals of YM2B2C (T c =15.5 K) were grown 
by the floating zone method. The residual resistivity ra- 
tio was approximately 47, indicating the highest crystal 
quality currently achievable. We measured the c-axis 
thermal conductivity k zz (the heat current q \\ c) on 
the sample cut into a rectangular shape (0.31 x 0.54 x 
3.00mm 3 ) by the steady-state method with one heater 
and two RuC>2 thermometers. To apply H with high ac- 
curacy relative to the crystal axes, we used a system with 
two superconducting magnets generating H in two mu- 
tually orthogonal directions and a 3 He cryostat equipped 
on a mechanical rotating stage with a minimum step of 
1 / 500 degree at the top of the Dewar Jl?], [ij] . Computer- 
controlling the two magnets and the rotating stage, we 
were able to rotate H with a misalignment of less than 
0.02° from the plane, which we confirmed by the simul- 
taneous measurement of the resistivity p zz . 

We first discuss the T- and H- dependence of k zz . 
Figure 1 (a) shows the T-dependence of k zz in zero 
field. Upon entering the superconducting state, k zz ex- 
hibits a tiny kink. The Wiedemann-Franz ratio L = 
i^zzPzz/T ~ 1.02Lo at T c is very close to the Lorenz 
number L = 2.44 x 10~ 8 ttW /K, indicating that the 



electronic contribution well dominates over the phonon 
contribution. The phonon peak at ~ T c /3 reported in 
LUNI2B2C was not observed , possibly due to the differ- 
ence in the phonon spectrum in the two compounds [B5| . 
The inset of Fig. 1 (a) shows the same data below 1 K. 
The T-dependence of k zz is close to quadratic rather than 
cubic. We will discuss this T-dependence later. Figure 
1 (b) depicts the .^-dependence of n zz (jy || [110]) at low 
temperaures. At low field k zz increases rapidly with H. 
This steep increase is markedly different from that ob- 
served in the ordinary s-wave superconductors in which 
the thermal conductivity shows an exponential behav- 
ior with much slower growth with H at H H C 2- It 
has been reported that in LuNi2B2C(T c =16.6 K), a close 
cousin of YNI2B2C, the independence of the in-plane 
thermal conductivity k xx is even steeper than UPt3 with 
line nodes jl2| . The steep increase of the thermal conduc- 
tivity in YM2B2C and LuNi2B2C, along with the \fll- 
dependence of the heat capacity 1 10 , O, are evidence 



that the thermal properties are governed by the delo- 
calized QPs arising from the nodes in the gap function. 
The theoretical understanding of the heat transport for 
superconductors with nodes has largely progressed dur- 
ing past few years p(|. The most remarkable effect on 
the thermal transport is the Doppler shift of the energy 
of a QP with momentum p (e(p) — ► e(p) — v s ■ p) in 
the circulating supercurrent flow v s j27|. This effect be- 
comes important at such positions where the local en- 
ergy gap becomes smaller than the Doppler shift term 
(A < v s ■ p), which can be realized in the case of super- 
conductors with nodes. This effect gives rise to a steep 
increase of the thermal conductivity with H, which is 
consistent with the present results. Having established 
the predominant contribution of the extended QPs in the 
thermal transport, the next question is the nodal struc- 
ture of the gap function. The important advantage of 
choosing to measure the thermal conductivity is that it 
is a directional probe, sensitive to the relative orienta- 
tion between the ma gnet ic field and nodal directions of 
the order parameter g ||, ||, || ||, || §| ||. 
This statement is based on the fact that the magnitude 
of the Doppler shift depends on the angle between the 
node and H . For instance, when H is rotated within 
the basal plane in <i-wave superconductors, the Doppler 
shift gives rise to the fourfold oscillation of the den- 
sity of states (DOS). In this case, the DOS shows the 
maximum (minimum) when H is applied to the antin- 
odal (nodal) directions. In fact, a clear fourfold mod- 
ulation of the thermal conductivity, which reflects the 
angular position of line nodes of d-wave symmetry, has 
been observed in high-T c cuprate YBa2Cu307_,s fT(i| . 
heavy fermion CeCoIn 5 |nj , organic k-(ET) 2 Cu(NCS) 2 
JbSt , demonstrating that the thermal conductivity can be 
a relevant probe of the superconducting gap structure. 



Figure 2 displays 
of n zz , which was 
H =H(s'm 9 cos 0, sin 9 sin </> 
of 4>, keeping 9 constant. 



the angular variation 
measured by rotating 
, cos 9) conically as a function 
Here 9=(q, H) is the polar 
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FIG. 2: Angular variation of the c-axis thermal conductivity 
k zz at H=l T and T=0.43 K (thermal current qr|| c). k zz are 
measured by rotating H—H (sin 6* cos 0, sin 6 sin 0, cos #) coni- 
cally as a function of <j> at fixed 6(— 90°, 60°, and 45°). 9—(q, 
H) is the polar angle and cj> is the azimuthal angle measured 
from the o-axis (see the inset). The crosses represent the 
data at 9 = 90° which are obtained under the field cooling 
condition at each angle. 



angle and <f> is the azimuthal angle measured from the 
a-axis (see the inset of Fig. 2). The measurements have 
been done by rotating (f> after field cooling at <j> = —45°. 
The crosses in Fig. 2 show k zz (H,4>) at H=l T which 
are obtained under the field cooling condition at each 
angle. k zz (H, <j>) obtained by different procedures of field 
cooling well coincide with each other, indicating that 
the field trapping effect related to the vortex pinning is 
negligibly small at H=l T. Below 0.5 T, a small field 
trapping effect was observed (see also Fig. 1(b)). A clear 
fourfold symmetry is observed at 9 = 90° and 60°; k zz 
can be decomposed as k zz = kP zz + K z f where k° zz is 
a (^-independent term and K^f is a term with fourfold 
symmetry with respect to 0-rotation. Similar results 
were obtained at 0.27 K and 0.8 K. We note that in the 
previous measurements in which the in-plane thermal 
conductivity k xx was measured in H rotating within the 
a6-plane, a large twofold term appears as a result of the 
difference of the effective DOS for QPs traveling parallel 
to the vortex and for those moving in the perpendicular 
direction [FJ, |l6[ 0, 
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In the present geometry, on 
the other hand, such a twofold term is absent because 
out-of-plane thermal conductivity k zz was measured in 
keeping 0, i.e. the angle between H and q, constant. 
This geometry enables us to make a precise analysis 
of the term with fourfold symmetry, which is directly 
related to the electronic structure |M . 

As seen in Fig. 2, K^f shows a peculiar angular varia- 
tion. The first point to emphasize is the appearance of a 
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FIG. 3: Angular variation of I zz {<j>,d) normalized by 
7 ZZ (45°,90°) as a function of <j> at 9 = 90° (solid lines), 60° 
(dashed lines), and 45° (dash-dotted lines), which are calcu- 
lated by the gap functions with (a) point node and (b)line 
node. For the definition of 9 and (j), see the inset of Fig 2. 
The corresponding gap functions are illustrated in the insets. 
For details, see the text. 



a narrow cusp at (f> — 0° and 90°. This angular variation 
is markedly different from those in e?-wave superconduc- 
tors with line nodes, in which the fourfold oscillation is 
close to a sinusoidal wave |l^, Hi], 17, [l8|. The second 
point is that the amplitude of refj decreases rapidly as H 
is changed from 8 = 90° to 45°. Here we stress that the 



fourfold anisotropics of the Fermi velocity vf and H C 2, 
which are inherent to the tetragonal band structure of 
YM2B2C, are quite unlikely to be an origin of the ob- 
served fourfold symmetry for the following reasons. We 
measured the ^-dependence of H C 2 at 6 — 90° and 45° 
and found that H C 2 shows nearly perfect sinusoidal (f>- 
dependence with fourfold symmetry, similar to H C 2 of 
LUN12B2C reported in Ref.|gq|. This is totally different 
from the (^-dependence of kjj displayed in Fig. 2. In 
addition, the amplitude of fourfold oscillation of H C 2 at 
9 = 45° was nearly 1/3 of that at 6 = 90°, while the 
amplitude of kjj at 6 = 45° is less than 1/5 of that at 
9 = 90°, indicating a different ^-dependence of two quan- 
tities. These results are consistent with the theory of the 
thermal conductivity. According to the calculation based 
on the Kubo formula, the leading term is the anisotropy 
of the gap function and anisotropy of vp will only enter 
as a secondary effect. Thus these consideration lead us to 
conclude that the observed fourfold symmetry originates 
from the superconducting gap nodes. 

We are now in position to discuss the nodal structure 
of YNi2B2C. The fact that K z f shows the minimum at 
= 0° and 90° immediately indicates that the nodes are 
located along [100] and [010] -directions. The cusp struc- 
ture and the ^-dependence of key features for 
specifying the type of nodes. In Figs. 3 (a) and (b), 
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we compare k 4 ^ f° r two different types of nodes calcu- 
lated from the Doppler-shifted QP spectrum. Here we 
adopted gap functions A(k) = hAn {l — sin 4 9 cos(4</>)} 
(s + g-wave) proposed in Ref. p2|, p3| for point node, and 
A(fe) = Ao sin(2</>) (d-wave) for the line node. These gap 
functions are illustrated in the insets of Figs. 3 (a) and 
(b). We caution readers that the the gap functions are 
rotated by 45° from those of Ref. (2^, |2^j, since the nodes 
in these references are located along the [110] direction. 
In H for general orientation given by (9, 0), k zz (H, <fi, 9) 
is given as 



where v ajC are the anisotropic Fermi velocities. For point 
nodes I zz is given as I zz (<fc, 9) = /+(</>, 9) 3 and 

I+((f>,9) = -[{1 + cos 2 6» + sin 2 6icos(20)}^ 

+{l + cos 2 6»-sin 2 6»cos(20)}5]. (2) 

For the line nodes I zz is given as I zz ((j),9) = I + (cj),9) 2 
and 

l + (<f>,0) = — / dV[{l + ^sin 2 6'(cos(2(/ ) )-cos(2^)) 

+ ^sin(20) simpy/l - sin(26»)}^ 
v2 

+{1 - - sin 2 6»(cos(20) + cos(2^)) 

+^sin(20) sin^v/1 + sin(20)}*].(3) 
v2 

Here, the superclean limit <C tP- is assumed, where 
r is the carrier scattering rate which is estimated to be 
HT <1 K [§§. Using A -28 K, # c2 ~10 T, the su- 
perclean condition is well satisfied in the present mea- 
surement. We first compare K^f at 9 = 90°. A clear 
narrow cusp structure is seen at <fi = 0° and 90° for the 
point node, while it is absent and K z f is close to the si- 
nusoidal wave for the line node. Qualitatively this can 
be explained as the following. The narrow cusp in K z f 
appears as a result of the smaller freedom of the QPs in- 
duced in the vicinity of point nodes. On the other hand, 
the cusp structure is smeared out due to the larger free- 
dom associated with the line node. We next discuss how 
the amplitude of changes with 9. For point nodes, 
the amplitude of K z f at 9 = 45° is much smaller than that 
at 9 = 90°, while they are of almost the same magnitude 
for line nodes. This can be accounted for considering the 
fact that rotating H conically at 9 = 45° does not point 
to the nodes in case of the point nodes, while rotating 
H at any 9 always cross the nodes in case of the line 
nodes. It is apparent that both the cusp structure and 
^-dependence of K^f shown in Fig. 3(a) are strongly in 
favor of the point nodes. We have calculated n z f for dif- 
ferent types of gap functions with point nodes and found 



that both the cusp structure and the ^-dependence of K z f 
are robust for the choice of the gap functions. In addi- 
tion, we have measured the in-plane thermal conductivity 
(q\\ [110] ) of a different crystal in H rotated within the 
[110]-planc and found no anomaly associated with the 
point node along the [001]-direction. These results lead 
us to our forementioned conclusion that the supercon- 
ducting gap function of YM2B2C has point nodes which 
are located along [100] and [010] -directions. 

For point nodes, the T-dependence of the thermody- 
namic quantities at low temperatures strongly depend 
on the dispersion of the gap function near the nodes. 
We point out that the gap function we used predicts 
a quadratic T-dependence of the thermal conductivity, 
which is consistent with the data in the inset of Fig. 1(a). 
On the other hand, a cubic T-dependence of k xx has been 
reported in LuNi 2 B 2 C below 150 mK pi pi. One pos- 



sible origin for the discrepancy between YN^B^C and 
LuNi2B 2 C is the difference of the dispersion of the gap. 

We here comment on the gap anisotropy and the im- 
purity effect. In the present experiment, we cannot rule 
out the possibility that there remains a very small but 
finite gap at the nodal positions. Let us estimate the up- 
per limit of the gap anisotropy ratio A m ,- n / A max , where 
A max and A m i n are the gap maximum and minimum, 
respectively. The fact that the clear fourfold pattern 
characterisitic to the point node was observed even at 
T=0.27 K suggests that A TO j„ is much less than 0.27 K. 
Assuming the weak coupling BCS relation for the maxi- 
mum gap (A max =l.75k B T c ), we estimate A max ~28 K. 
We then obtain A m j„/ A max <C 0.01; the anisotropy ratio 
is essentially zero. According to Ref. [[l0|, E3, the intro- 
duction of non-magnetic impurity in YM2B2C reduces 
the DOS at the Fermi surface by removing the node. 
This indicates that the gap function does not change its 
sign at the nodes, because otherwise the impurity always 
induces a finite DOS at the Fermi level, as observed in 
high T c cuprates |3(J. Therefore the impurity effect is 
consistent with the point node for which case the gap 
function does not change its sign at the node. 

Various authors have pointed out that calculations 
based on the Eliashberg equation reproduces well the 
high transition temperature of borocarbide supercon- 
ductors in terms of the observed phonon spectrum and 
isotropic gap |J. However the determined gap struc- 
ture with point nodes casts doubt on such view based on 
a simple phonon mediated pairing mechanism. At the 
present stage, what kind of the pairing interaction gives 
rise to the point node is not known. We also point out 
that the determined nodal structure should play an im- 
portant role in determining the superconducting proper- 
ties, such as the vortex lattice structure, reversible mag- 
netization, upper critical field H C 2, etc. Until now, these 
preperties have been discussed in terms of an ordinary s- 
wave superconductor with anisotropic Fermi surface with 
the use of nonlocal London theory p3| . However, the 
determined nodal structure motivates further investiga- 
tions on these subjects. For instance, the extended QPs 
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appear to be very important for the vortex triangular- 
square lattice phase transition []3l"|| . 

In summary, we demonstrated that the measurement 
of k zz in H rotating conically around the c-axis with a 
fixed tilt angle provides a novel method to distinguish line 
and point nodes. We presented the first compelling evi- 
dence that the gap function of YM2B2C has point nodes 
which are located along the [100] and [010] directions. To 
the best of our knowledge, YN^E^C ( and presumably 



LUM2B2C) is the first superconductor in which the gap 
function with the point nodes was successufully identi- 
fied. The determined gaps structure offers a new per- 
spective on the pairing mechanism as well as unusual 
superconducting properties of anisotropic superconduc- 
tors. 
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